Shapiro steps in Josephson junctions with alternating critical current density 
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We treat theoretically Shapiro steps in tunnel Josephson junctions with spatially alternating crit- 
ical current density. Explicit analytical formulas for the width of the first integer (normal) and 
half-integer (anomalous) Shapiro steps are derived for short junctions. We develop coarse-graining 
approach, which describes Shapiro steps in the voltage-current curves of the asymmetric grain bound- 
aries in YBCO thin films and different superconductor-ferromagnet-superconductor Josephson-type 
heterostructures. 
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I. INTRODUCTION 

Series of resonances exist in Josephson tunnel junc- 
tions biased at an alternating voltage V{t) — + 
Vi cos wii.— "^iS In junctions of conventional superconduc- 
tors these resonances appear at the frequencies loq — 
2eVo/h — nuji, where n is an integer. At the resonant 
values of the dc voltage Vb — nfkui /2e the supercurrent 
has a dc component. In the voltage-current curve Shapiro 
resonances reveal itself as a "ladder" of equidistant values 
of Vo {integer Shapiro steps) fii^»2. 

The physical origin of Shapiro steps follows from the 
Josephson equations j = jc sin ip and ip = 2eV/fi, where 
j and jc are the tunneling and critical current densi- 
ties, and (f is the phase difference across the junction. 
In order to find the current across the junction we inte- 
grate the voltage V{t) and find the time-dependent phase 
(f{t) — (^0 + LO{)t + vi sinciJii. In this relation ipa = (/3(0) 
is the initial value of the phase and vi = 2eVi /fiLOi is 
the dimensionless parameter of the problem. Knowing 
Lp{t) we obtain the tunneling current density in the form 
j(t) = jc sin((/Jo + ojQt + vi siwujit). This formula demon- 
strates that j{t) is a complex alternating function of time. 
Fortunately, j{t) can be transferred into a series allowing 
for an easy qualitative and quantitative analysis^i^ 
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(-1)" J„(wi) sin[(/7o + (wo - nuJi)t], (1) 



where Jn{x) is the first kind Bessel function of order n. 

It is seen from Eq. ([1]) that for any resonant frequency 
Wo = noji the supercurrent density j has a dc component 
oc sin(/?o, which reaches its maximum, jm, at (^o = 7r/2. 
As a result the maximum value of the dc current density 
is given by jm = jcJn{vi) ^ 

Anomalous Shapiro steps at the subharmonic reso- 
nant frequencies wq = {n/q)LOi (q > n, where q and n 
are integers) were treated theoretically for short micro- 
bridges of conventional superconductors assuming that 
the current-phase relation includes high-order harmonic 
terms.* Qualitatively the effect of terms jq sm{qip) on the 
supercurrent density j follows from Eq. ([1]) . Indeed, sub- 
stituting ifQ, ujQ and vi by q^po, qtoo and qvi we find that 
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FIG. 1: Schematic diagrams: (a) Josephson junction made 
of periodically interchanging - and tt -biased fragments; (b) 
alternating critical current density distribution. 



at any subharmonic frequency ljq = {n/q)u>i the super- 
current density j has a dc component oc svaqtpQ. In this 
case the maximum value of j corresponds to the phase 
V?o = 7i'/2(j and is given by jm = jqJn{qvi). 

Observation of anomalous half-integer Shapiro steps 
has been reported recently for asymmetric grain bound- 
aries in YBCO films^ and superconductor-ferromagnet- 
superconductor (SFS) heterostructures.i^i^ 

Qualitative explanations of the origin of the anomalous 
Shapiro steps were proposed to understand the results of 
these experiments mainly by revealing the existence of 
second harmonic term, sin 2ip, in the current-phase re- 
lation (CPR). It was shown that data measured for the 
asymmetric grain boundaries in YBCO are indeed consis- 
tent with existence of the term sin 2ip in the CPR^ '^°i^^ 
In the case of SFS junctions, it was assumed that there 
is doubling of the Josephson frequency, which leads to 
a sin2(y9 term in the CPR. It was suggested that this 
frequency doubling is caused by splitting of energy lev- 
els in the ferromagnetic exchange fieldi^ An alternative 
explanation^ assumes existence of a resonance between 
the spontaneous currents and rf modulation. 

Asymmetric grain boundary junctions in YBCO thin 
films and SFS heterostructures are arranged in sequences 
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of interchanging 0- and tt -shift tunneUng junctions as 
shown in Fig. [TJ tt -shifts spontaneously appearing 
across superconducting banks of tunnel junctions were 
first considered for SFS heterostructuresJ^ii^ These tt- 
shifts lead to negatively biased critical current density, 
jc, and anomalous Josephson properties of sequences 
of O-TT biased tunnel junctionsJ^ii^iilii^iiaiSa The ob- 
served anomalies are especially significant if the lengths 
of - and tt -fragments are much less than the Josephson 
penetration depth Aj defined by the spatial average of 

|„- I 21 .22.23.24 
\jc\' 

Thus, several qualitative explanations have been put 
forward to understand the physics of fractional Shapiro 
steps. However, to the best of our knowledge there is 
at present no theoretical model describing the anoma- 
lous Shapiro steps and their basic features in Josephson 
junctions with alternating critical current density. 

In this paper we treat Shapiro steps in Josephson junc- 
tions arranged in periodic or almost periodic sequences 
of interchanging 0- and tt -biased fragments with I <ti Xj. 
We develop a coarse-graining approach which is applica- 
ble to treat time-dependent phenomena. 

In our theoretical model we assume that the tunnel- 
ing current density is given by the standard CPR depen- 
dence j = jcix) smip{x,t), but with a critical current 
density, jc{x), spatially alternating along the Josephson 
junction. As an illustration we treat the case of a short 
junction that allows to obtain explicit analytical results. 
We demonstrate that in Josephson junctions with spa- 
tially alternating critical current density the half-integer 
(anomalous) Shapiro steps exist in addition to the integer 
(normal) Shapiro steps. We calculate the dependence of 
the width of these steps on the flux inside the junction 
and the voltage across the junction. 

The organization of this paper is as follows. In Sec. II 
we review briefly the coarse-graining approach to Joseph- 
son junctions with spatially alternating critical density. 
In Sec. Ill this general approach is applied to calculate ex- 
plicit formulae for the width of the integer (normal) and 
half-integer (anomalous) Shapiro steps in short junctions. 
Sec. IV discusses these results. 



II. JOSEPHSON JUNCTION WITH SPATIALLY 
ALTERNATING CRITICAL CURRENT DENSITY 

Consider a tunnel junction of length L {0 < x < L) 
with a critical current density jc{x) alternating with a 
length-scale / ^ A, where A is the London penetration 
depth (see Fig. [T]for the geometry of the problem). It is 
convenient for the following calculations to write 

Jc{x)^{m+9{x)], (2) 

where (jc) is the average value of jc(x) defined as 



(/) = - J^dxfix). 



(3) 



The function g{x) describes spatial variations of jc{x) 
and has a zero average value, {g{x)) = 0. In what follows 
we are interested in the case of small average values of the 
critical current density, i.e., we assume that {jc{x)) 
{\jcix)\). This condition means that the typical value of 
the dimensionlcss function g{x) is big compared to unity 
or in other words (|g(x)|) ^ 1. 

In the above notations equation for the phase difference 
(p{x, t) across the junction reads 



AV" + [^ + 9{x)] sin LP ^0, 



(4) 



where l/r is the Josephson frequency, and A is the effec- 
tive Josephson penetration depth, 



A^ = 



C0O 



(5) 



We treat now a sample subjected to magnetic field Ha 
and alternating voltage V{t) = Vq -I- Vicoswii applied 
across the junctions. It is assumed that 2eVi is small 
compared to tkoi, i.e., vi = 2eVi /TiliJi <C 1. In this case 
coarse-graining can be applied to solve Eq. (|4]).5S 

Two types of terms appear in Eq. ([4]): terms alter- 
nating over the length / and smooth terms varying over 
the length A ^ /. The fast alternating terms cancel each 
other, independently of the smooth terms, which also 
cancel each other. Therefore, we write^^ 



(p{x) = ^p{x) +^ix). 



(6) 



where ipix) is a smooth function with the length-scale of 
order A and £,{x) alternates with the length-scale of order 
I. Under the above assumptions the average value of ^{x) 
is zero and the typical amplitude of variations of ^(x) is 
smah. I.e., (C(a;)) = and {\({x)\) < 

Substituting Eq. ^ into Eq. ^ and keeping terms up 
to first order in ^(a;) we findS^ 



A2^" - 

(jc) 

2e// jd^) 



= 0, 



0, 



(7) 

(8) 



where the smooth j^{x) and alternating j^{x) compo- 
nents of the tunneling current density j — + are 

ji: = (jc) (sin?A - 7sinV'C0sV') , (9) 

j? = Uc) 9{x)smil} . (10) 

The dimensionlcss constant 7 is equal to 

7 = (5(^)G(^)>-A2(e;2(x))>0, (11) 

and the rapidly alternating phase ^g{x) is defined by 

ax) = -U^)sin4,. (12) 

It follows from Eqs. ©, (dH]) and HID that 



AX+9{x) = 0, 



(13) 
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i.e., the rapidly alternating phase shift ^g{x) is a charac- 
teristics of a sample. 

To summarize the coarse-graining approach it is worth 
noting that the typical values of S,g{x) are small, but at 
the same time the typical values of g{x) are big, i.e., 
{\^gi^)\) 1 (|g(a;)|) 1. As a result, the dimen- 
sionless parameter 7, which is proportional to the average 
of the product of the two rapidly alternating functions 
£,g{x) and 17(2;) might be of the order of unityi^i^ The 
value of can be estimated as ~ Cg/^- It follows then 
from Eq. ^ that (^^^ ~ jP/A^ < 1. Therefore for 
we have the estimate A(|^^|) ~ ^ ^ 1. 

III. SHAPIRO STEPS IN SHORT JUNCTIONS 

In a short junction {L <^ A) 'il'{x, t) is almost linear in 
X, the time dependence of ^{x, t) is given by V{t), and 



00 

i^{L/2,t)^'^-^V{t), 



(14) 
(15) 



where Hi(t) is the time-dependent field in the junction. 

At the sample edges the derivative ip' = ip' + ^' is pro- 
portional to the field i?|o,L including the self-field gen- 
erated by the total current /jSi In the case of a short 
junction we obtain 

w^It) -viox-Avy|„,,+A^'l„^, (16) 

where Hg = (po/inXA, Is = (jc)A. It follows from Eq. 
(fT^ that the derivative ^' = — sin V'. In order to find 
we integrate Eq. (|13p from to L and arrive to the rela- 
tion Cgi^) = Cg{L), which allows to rewrite the boundary 
condition in the final form 



Ha 

Hs 



1 / 
2Z 



= Ai/i'lg L — a sini/j 



0,L 



0,L ■ 



(17) 



where the terms a sin "0 1 q l caused by the high-density 
edge currents and a is a constant. It is worth noting here 
that a is a characteristics of a sample. Next, using the 
above estimate for we obtain 



(18) 



To summarize the above analysis, we find that in the 
coarse-graining approach the Josephson junctions with 
spatially alternating critical current density and / ^ Aj 
are characterized by two dimensionless parameters 7 (for 
the inner part of the junction) and a (for the edges of 
the junction), where Xj = ^yc(j)o/16^T'^X{\jc\) <C A is the 
local Josephson length. 

Next, we find the derivatives '/''Ig l ^^i^S the first in- 
tegral of Eq. ([7]) , that can be written as 



A' 



2 V'^ + cos^l!-J cos 2-0 = Const. (19) 



Finally, we combine Eqs. ([T7|) and (fTQ]) relating the total 
current /, applied field Ha, and phase at the edges. In 
the case of a short junction we find 



I = 2Is sin I TT 
L 



27rA 



1 

a L 



a cos 1pm - 
f 

7 cos 0m COS TT ■ 



siniPm M20) 



-0" 



'0 , 

sin -0m, (21) 
ip{L/2,t) = tjjo + ujot + visinojit, (22) 



where the "applied", 0q, and "internal", 0^, fluxes are 
defined as 0a, i — 2LXHa.i, ujq = 2eVo/h, and ipQ is a 
constant, which is used to maximize the total current /. 

Calculation of the width of Shapiro steps similar to the 
one given by M. Tinkham^ leads to the following results: 
two series of steps appear at frequencies 



ujQ — (n + 1/2) LOi and (jJo = nuii, 



(23) 



where n is an integer. 

In the case of low oscillating voltage {vi ^ 1) we obtain 



the widths of the first half-integer, Ii {(pa 
integer, /i(0a), Shapiro steps in the form 

700 



and the first 



Ii = Vilc 



2 2 
a cos 



2tt(i 



sin I 27r— 

00 




(24) 



(25) 



where Ic = {jc)L. Explicit formulas ([24]) and ([25]) reveal 
quite a few remarkable features of Ikiipa) and /i(0a). In 
what follows we discuss them in details. 



IV. RESULTS AND DISCUSSION 

The function Ii{<j>a) is equal to zero for the two se- 
ries: 0a = (« + 1/2) 00 and 0a = (n — ■06//7r)0o, where n 
is an integer and ipe is defined by tan ipg = ira^cpa/jft'o- 
The angle ipe depends on 0a, i.e., in general, the roots of 
equation Ii (0a) = are not equidistant and the actually 
observed dependence /i (0a) is defined by the dimension- 
less ratio ipa = 7/27ra^. 

Next, it follows from Eq. that for 0a = the width 
of the first anomalous half-integer step Ii (0) equals to 
vilc{l + o?). In the low- flux region (0a /0o ^ V'a) the 
main contribution to Ii{ipa) comes from the second term 
in Eq. (|24p . This term originates from the alternating 
currents fiowing across all junctions. As a result, in the 
low-flux region we have Ii oc (0o/27r0a)| sin(27r0a/0o)|, 
i.e., /i(0a) is described by the Fraunhofer pattern but 
with a double frequency. In Fig. [l](a) we show /i (0a) in 



4 




FIG. 2: Dependence of the width of the first half-integer 
Shapiro step, Ii, on applied fiux, (j)a, given by Eq. (|24p for 
7 = 1.5 and (a) a = 0.03, (b) a = 0.3, and (c) a = 1. 



the low-flux region {4>a/4'o 250) for the data a = 0.03 
and 7 = 1.5, which lead to ipa ~ 250. 

In the high-flux region {(pa/ipo ^ V'a) the main contri- 
bution to /i comes from the first term in Eq. ([24| . This 
term originates from the high density currents flowing 
across the edge junctions. As a result, in the high-flux 
region we have /i oc cos'^{TTipa/4>o), i.e., I^{(f>a) is a pe- 
riodic function. In Fig. [2](c) we show /i(0a) in the 
high-flux region (0a/0o ^ 0.25) for the data a = 1 and 
7 = 1.5, which lead to ipa ~ 0.25. 

In Fig. [11(b) we show the function Ii {(j)a) for the inter- 
mediate values of the applied flux {(j)a/4'o V'a)- Using 
the data 7 = 1.5 and a = 0.3, we find that V'a ~ 2.5. In 
this case ij^g is strongly flux dependent. As a result the 
fimction Ii{(l)a/4'o) is manifestly aperiodic. 

We discuss now the width of the first integer Shapiro 
step, /i. It follows from Eq. that /i((/'a) is equal to 
zero for the series: 4'a = rupQ, where n is an integer. 

In the high- flux region [(j)a/4>o 3> Lji^aK) the main 
contribution to 1\ comes from the first term in Eq. (|25p . 
This term originates from the high density currents flow- 
ing across the edf^e junctions. As a result, in the high-flux 
region we obtain l\{<^a) oc | sin(7r0a/</'o)|, i-e., /i(0a) is 
a periodic function. In Fig. [31(a) we show /i((/'a) in the 
high- flux region {4)al4>o ^ 0.2) for the data a = 0.3 and 
L/A = 0.2, which lead to L/iraK w 0.2. 

Next, it follows from Eq. (|25p that the width of the 
first integer step /i(0) equals to vilc- In the low- flux 
region (jpa/ipo ^ L/iraA) the main contribution to Ii 
comes from the second term in Eq. (j25p . This term 
originates from the alternating currents flowing across all 
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0.5- 



FIG. 3: Dependence of the width of the first integer Shapiro 
step, 7i, on applied flux, 0a, given by Eq. (|25p for L/A = 0.2, 
and (a) a = 0.3, and (b) a = 0.003. 



junctions. As a result, in the low-flux region we obtain 
Ii oc ((/)o/7r(?!)a) sin(7r0a/0o), i.e., Iii4>a) is described by 
the Fraunhofer pattern. In Fig. [21(b) we show Ii{(l>a) in 
the low-flux region {<j)a/<j)Q ^ 20) for the data a — 0.003 
and L/A = 0.2, which lead to L/iraA « 20. 

Let us now illustrate the above calculations by using a 
model dependence for alternating critical current density 



■ , ^ ■ . . / 27r 

Jc{x) = Jo +Jism 1 — 



(26) 



where jo = {jc{x)) and ji are constants, 9 is an angle 
from the interval < 9 < tt, and L/l — N is an integer 
(iV > 1). It follows from Eqs. Q and ^ that 



Ji . /27r 
■ sm I ~ 2; 



Jo 
c4>o 



(27) 
(28) 



Next, knowing g{x) and ^g{x) we use Eqs. (fTT|) and (|20l) 
and find the parameters 



c4>ojo ' 
■\/27 COS 9 . 



(29) 
(30) 



It follows from Eqs. ([M]) and ([25]) that the widths of 
the first half-integer Ii and integer Ii Shapiro steps are 



vijiL j|_ 
cosipe Af 



cos 



COS 



jih 
An 



~L '^t 

00 



00 



(31) 



, (32) 



where we define tan^/^g = {(f)f)/2Tr(t)a)/ cos^ 9, li — lcos9, 
and Al = c(/)o/4Aji. 
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An interesting feature follows from Eq. ([32]) for the lim- 
iting case jo = 0- Using Eg. ([5^ we find that if the aver- 
age value of the critical current density jo ~ 0, then the 
width of the first integer Shapiro step equals to 



— — — cos 



sm 



4>o 



(33) 



This current is the contribution coming from the edge 
fragments and therefore the value of Ii does not depend 
on the total length of the junction L. 



V. SUMMARY 

The anomalous Shapiro steps in our model exist due 
to: (a) successful interference between the spatial alter- 
nations of the critical current density jc{x) and phase 
factor sin f{x, t), leading to generation of the second har- 
monic in the Josephson current density; (b) high current 
density at the edges of the junction resulting in anoma- 
lous dependencies Ii{Ha) in terms of periodicity at low 



fields and asymptotic behavior at high fields. 

To summarize, we demonstrate the existence of anoma- 
lous half-integer Shapiro steps for short Josephson junc- 
tions with spatially alternating critical current density 
ic{x). We derive exphcit formulas given by Eqs. (|24|) 
and ([25|) for the width of the first integer and half-integer 
Shapiro steps. The general approach is applied to the 
case of a simple model for the critical current density 
dependence on the coordinate along the junction (see 
Eq. (Uni))- The results obtained in the framework of 
this model might be useful for analysis of the experimen- 
tal data for the asymmetric grain boundaries in YBCO 
thin films and different superconductor-ferromagnet- 
superconductor Josephson-type heterostructures. 
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